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Vortex trapping and expulsion in thin- film YBa2Cu307_5 strips 
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A scanning SQUID microscope was used to image vortex trapping as a function of the magnetic 
induction during cooling in thin-film YBa2Cu307_5 (YBCO) strips for strip widths W from 2 to 50 
/im. We found that vortices were excluded from the strips when the induction Ba was below a critical 
induction Be. We present a simple model for the vortex exclusion process which takes into account 
the vortex - antivortex pair production energy as well as the vortex Meissner and self-energies. This 
model predicts that the real density n of trapped vortices is given by n = [Ba ~ Bk)/^o with 
Bk = 1.65^0 /W^ and $o = h/2e the superconducting flux quantum. This prediction is in good 
agreement with our experiments on YBCO, as well as with previous experiments on thin-film strips 
of niobium. We also report on the positions of the trapped vortices. We found that at low densities 
the vortices were trapped in a single row near the centers of the strips, with the relative intervortex 
spacing distribution width decreasing as the vortex density increased, a sign of longitudinal ordering. 
The critical induction for two rows forming in the 35 fim wide strip was (2.89 -I- 1.91 — 0.93)-Bc, 
consistent with a numerical prediction. 

PACS numbers: 74.25.Ha, 74.25.Qt, 74.25.0p, 74.78. Bz 
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I. INTRODUCTION 



In principle, when a parallel magnetic field is applied 
to an infinitely long, defect-free superconducting cylin- 
der, all magnetic flux should be expelled as the tempera- 
ture T is lowered through the superconducting transition 
temperature Tc, provided that the applied magnetic field 
is below either the critical field Hc{T) for a type-I super- 
conductor, or the lower critical field Hci{T) for a type-II 
superconductorji In practice, real samples have finite size 
and often contain defects, which can pin magnetic flux. 
Moreover, nonellipsoidal samples, even those not contain- 
ing defects, naturally possess geometric energy barriers 
that can trap magnetic flux during the cooling process. 
Pinned or trapped vortices are nearly always observed 
in thin-film type-II superconductors, even when cooled 
in relatively low magnetic fields. In general, this can be 
attributed both to pinning of vortices by, for example, 
defects and grain boundaries, and to trapping by the geo- 
metric energy barriers. Understanding such pinning and 
trapping effects is important for superconducting elec- 
tronics applications. 

The present work is motivated by applications of high- 
Tc superconducting sensors such as SQUIDs^ and hybrid 
magnetometers based on high-Tc flux concentrators.'^ 
These sensors are used in a broad field of applications, 
such as geophysical research^ and biomagnetismi^ The 
sensitivity of these sensors is limited by 1// noise in an 
unshielded environment. The dominant source of this 
noise is the movement of vortices trapped in the sensor. 
This noise can be eliminated by dividing the high-Tc body 



into thin stripsi2i£ The strips have a certain critical in- 
duction below which no vortex trapping occurs, resulting 
in an ambient field range in which these sensors can be 
effectively operated. We investigated vortex trapping in 
thin-film YBCO strips in order to incorporate the results 
in a hybrid magnetometer based on a YBCO ring tightly 
coupled to, for example, a GMR (giant magneto resis- 
tance) or Hall sensor. 

Models for the critical induction of thin-film strips have 
been proposed by Cleini and Likharevji^. Indirect exper- 
imental testing of these models was done by observing 
noise in high-Tc SQUIDs as a function of strip width 
and induction.^'^'^ The induction mentioned here is the 
magnetic induction during cooling, which is the notation 
throughout this paper. More direct experimental veri- 
fication of these models was presented by Stan et al.— 
using scanning Hall probe microscopy (SHPM) on Nb 
strips. Both experiment and theory found that the crit- 
ical induction varied roughly like 1/W'^. However, the 
experimental^*^ and theoretical^i^ pre-factors multiplying 
this dependence differed significantly. In this pa- 

per we propose a model for vortex trapping in narrow 
superconducting strips which takes into account the role 
of thermally generated vortex-antivortex pairs. 

To test this model we performed scanning SQUID 
microscopy (SSM)l^ measurements on thin-film YBCO 
strips. We found excellent agreement between the depen- 
dence of critical induction on strip width and the present 
model for both our experiments on YBCO and for the 
previous work on Nb. In agreement with this previous 
work and as predicted by the present model, we found 
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that in YBCO the number of vortices increased for in- 
ductions above the critical induction linearly with the 
difference between the applied induction and the criti- 
cal induction. In a follow-up to the paper of Stan et 
al., Bronson et ali^ presented numerical simulations for 
the vortex distribution in narrow strips. These simula- 
tions showed that for inductions just above the critical 
induction the vortices are trapped in the centers of the 
strips. For higher inductions the vortices formed more 
complex ordered patterns, first in two parallel rows, then 
for higher inductions in larger numbers of parallel rows. 
We performed statistical analysis of the vortex distribu- 
tion in our measurements and found agreement with this 
model. 



II. THEORY OF VORTEX TRAPPING IN A 
THIN FILM STRIP 

Whether or not a vortex gets trapped in a strip is de- 
termined by the Gibbs free energy. This energy exhibits 
a dip in the center of a superconducting strip for applied 
inductions above a certain critical value. This dip gives 
rise to an energy barrier for the escape of the vortex. The 
models proposed by Clem^ and Likharev^ differ from the 
present model only in the minimum height of the energy 
barrier required to trap vortices. 



A. The Gibbs free energy of a vortex in a strip 

Consider a long, narrow, and thin superconducting 
strip of width W in an applied magnetic induction Ba- 
The vortex trapping process occurs sufficiently close 
to the superconducting transition temperature that the 
Pearl length A — 2X^ /d, with A the London penetration 
depth and d the film thickness, is larger than W. In 
this limit there is little shielding of an externally applied 
magnetic induction Ba- The resultant superconducting 
currents in the strip can be calculated using the fluxoid 
quantization condition: 



in response to the applied induction is: 



B ■dS + fioX^ (h Js-ds^ N<i>o. 



(1) 



In this equation the first integration is over a closed sur- 
face S within the superconductor, the second is over a 
closed contour surrounding S", _B is the magnetic induc- 
tion, Jg is the supercurrent density, and N is an inte- 
ger. SI units are used throughout this paper. If we take 
the strip with its long dimension in the y direction, with 
edges at X = and x — W, and an applied induction 
perpendicular to the strip in the z direction, a square 
closed contour can be drawn with sides at y — ±1/2 and 
X — W/2±Ax. If we assume uniform densities n-^ and 
of vortices and antivortices in the film, with n = — Ua 
being the excess density of vortices over antivortices, the 
first integral in Eq. ^ becomes 2BaAxl, the second be- 
comes 2Jsl, N = 2nlAx, and the supercurrent induced 
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{Ba - n$o) [x - WI2) . 



(2) 



The assumption of a uniform density of vortices is good 
at high trapping densities, and at zero density, but is in- 
correct at low densities, as we shall discuss later. Equa- 
tion ([2]) differs from the expression given in Ref. by the 
term — ri^o- As vortices are nucleated in the film, they 
reduce the screening currents induced in the filrn by the 
applied induction. The equation derived in Ref. {J\ for the 
Gibbs free energy of an isolated vortex (upper sign) or 
an antivortex (lower sign) at a position x inside the strip 
is then slightly modified as :^^'^^ 
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The Gibbs free energy consists of two terms. The first 
term, which is independent of the applied magnetic in- 
duction Ba, is calculated to logarithmic accuracy, as it 
includes only the kinetic energy of the supercurrents, and 
it is equal to $o-^circ/2, where Idrc is the supercurrent 
circulating around the vortex. This term, which has a 
dome shape and decreases monotonically to zero as the 
vortex reaches a distance ^/2 from the edges of the strip, 
is also equal to the work that must be done to move the 
vortex from its initial position ata; = ^/2ora; = W — ^/2 
to its final position at x against the Lorentz forces of at- 
traction between the vortex and an infinite set of neg- 
ative image vortices at —x + 2mW, m — 0, ±1,±2,... 
. Here (, is the coherence length, which is assumed to 
obey W. We also assume that the vortex core ra- 
dius is ^, such that the constant a = 2/7r as in Ref. 0. 
Other values of a, such as I/tt as in Ref. [l^, or 1/4 as in 
Ref.fi correspond to different assumptions regarding the 
core size. The second term in Eq. ([3|) is the interaction 
energy between a vortex (upper sign) [or an antivortex 
(lower sign)] and the screening currents induced by the 
external magnetic induction. It is the negative of the 
work required to bring a vortex (or antivortex) in from 
the edge against the Lorentz force due to the induced 
supercurrent given in Eq. ([2]) . The upper sign in Eq. Q 
corresponds to the fact that Jay tends to drive vortices 
into the film, and the lower sign indicates that antivor- 
tices are driven out. When Ba is sufficiently large, this 
term makes a minimum in G{x) in which vortices can be 
trapped. For wider strips this minimum occurs at lower 
values of the induction. 



B. Previous models for the critical induction 

There are two existing models which predict the crit- 
ical induction for vortex trapping when the applied per- 
pendicular magnetic induction is small {Ba ~ ^o/W^). 
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In these models the Gibbs free energy from Eq. ([3]) is 
used in the Umit of n — > 0. The critical induction model 
by Likhare"\*^ states that in order to trap a vortex in a 
strip the vortex should be absolutely stable. This hap- 
pens when the Gibbs free energy in the middle of the 
strips equals zero and leads to 



a vortex-antivortex pair is given by the pairing energyji 
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where a is the constant in Eq. 

Another model for the critical induction is proposed 
by Clem^i who considers a metastable condition. In 
this view vortex trapping will occur when the applied 
magnetic induction is just large enough to cause a mini- 
mum in the Gibbs free energy at the center of the strip, 
(fG{W/2)/dx^ = 0, leading to 
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the vortex-antivortex pair generation rate is given by a 
pre-factor times the Arrhenius factor exp{— Epair / ksT) , 
where ks is the Boltzmann constant. The vortex escape 
rate is given by an attempt frequency times a second 
Arrhenius factor exp(— E's/fc^T), where Eb is the differ- 
ence in the Gibbs free energy between the local maximum 
and the minimum in the center of the strip. Since Eb 
and Epair have the same temperature dependences (recall 
that 1/A is proportional to Tc ~ T), the vortex genera- 
tion rate and its rate of escape will be exactly balanced 
at all temperatures (aside from a logarithmic factor in 
the ratio of the two pre-factors) when Eb and Epair are 
equal. This occurs at a critical magnetic induction Bk 
which is the solution of the equation 
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C. Our model for the critical induction 

The model proposed here is intermediate between the 
models presented in Ref. 't' and 's'. As the strip is cooled 
just below the superconducting transition temperature 
Tc, thermal fluctuations cause the generation of a high 
density of vortex-antivortex pairs. Similar to the pro- 
cesses determining the equilibrium densities of electrons 
and holes in semiconductors, the equilibrium densities 
of vortices and antivortices very near Tc are determined 
by a balance between the rate of generation of vortex- 
antivortex pairs, the rate of their recombination, and 
the rates with which vortices are driven inwards and an- 



which leads to the condition 
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where the maximum value of the left-hand side of the 
equation is taken with respect to x. This equation can 
be solved numerically, resulting in 
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It appears to be an interesting numerical coincidence that 
the solution to this equation gives a pre-factor (1.6525) 



tivortices are driven outwards by the current Jy. Ac- that is only different from tt /6 by about 0.5%. 



cordingly, very close to Tc, the densities and Ua of 
vortices and antivortices equilibrate such that their dif- 
ference n = n„ — Tia is very nearly equal to Ba/ (j^o and 
the current Jy [see Eq. ([2])] is practically zero. When 
Bq < Ba < Bl, \i \s energetically unfavorable for vor- 
tices and antivortices to be present in the strip, and as 
the temperature decreases and the energy scales of the 
terms in Eq. ^ increase, the densities of both vortices 
and antivortices decrease. While vortices and antivor- 
tices continue to be thermally generated, the antivortices 
are quickly driven out of the strip by the combination 
of the self-energy and field-interaction energy [note the 
lower sign in Eq. ([3])]. The antivortex density thus be- 
comes much smaller than the vortex density and becomes 
so small that the recombination rate is negligible. The 
value of n « n„ drops below Ba/4>o. Although when 
Bo < Ba < Bl it is energetically unfavorable for a vor- 
tex to be present in the strip, the vortex's Gibbs free 
energy has a local minimum at the center of the strip 
and the vortex must overcome the energy barrier before 
it can leave the strip. Since the energy required to form 



As the temperature decreases, A decreases and be- 
comes much less than W . This means that once the 
vortices are trapped in the local minimum and arc clus- 
tered around the middle of the strip, the potential well in 
which they sit changes shape. Recall that the calculation 
of both terms in Eq. ([3]) assume that A is larger than W . 
It would not even be possible to magnetically image the 
vortices in the vicinity of the freeze-in temperature be- 
cause the local field produced by each vortex is then so 
spread out. However, as the temperature decreases, the 
number of trapped vortices per unit length remains fixed 
and the applied field remains constant. For n >> 
the distribution of vortices (averaging over the intervor- 
tex spacing) takes on a dome-like shape, and vortex-free 
zones appear at the edges of the strip. The z component 
of the local magnetic induction should then be described 
by the equations given in Sec. 2.1 of Ref. [H. The field 
distribution in a single strip containing a central vortex 
dome in which the current density is zero is closely re- 
lated to the field distribution of a pair of parallel coplanar 
strips with a gap between them .^^i^^ 
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D. Behavior above the critical induction 



Because in the present model the screening-current 
density [Eq. ([2])] and the Gibbs free energy [Eq. ^] de- 
pend on n, the areal density of vortices (when no an- 
tivortices are present), we can expect that for apphed 
inductions Ba weh above the critical induction Bk the 
balance between the rates of vortex generation and es- 
cape occurs when 



Ba-n^n = BK = 1-65-^, 



(10) 



which can be inverted to give the density n of trapped 
vortices as a function of applied induction, 



Ba — Bk 



(11) 
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However, for Ba just above B^:, where n <C Bx/^o, 
we need to take into account the interactions between 
vortices more carefully. As shown by Kogan^i^ the inter- 
action energy between a vortex at (x, 0) and another at 
{xi, tji) in a strip of width W > h. — /d and thickness 
d is 
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(12) 



where x = t:x/W , Xi = irxi/W , and iji = nyi/W. To ob- 
tain the interaction energy of the vortex at {x, 0) with 
all the vortices when the vortex density is very large 
(n <C l/W"^), it makes sense to convert the sum over all 
i to an integral, assuming a uniform density n over the 
strip width. The integral can be evaluated analytically, 
and the result is exactly equal to the term associated 
with n$o in Eq. ([3]). However, when the average vortex 
density is small, we can obtain an approximation to the 
interaction energy of the vortex at (a:, 0) with all other 
vortices by again converting the sum to an integral, as- 
suming uniform density of vortices n over the strip width 
but only for \yi\ > \/2nW. In other words, we exclude 
from the integral a rectangular region of height h and 
width W around the origin associated with one vortex, 
where n = 1/Wh. After changing variables of integra- 
tion, the resulting interaction energy Uint{W/2) at the 
center of the strip can be expressed as 
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where 0„ — sin ^[1/ cosh(7r/2nW^)]. For a given vortex 
density n, the interaction energy Uint{W/2) is less than 
the value that would be expected from the second term 
in Eq. We may define an effective vortex density 

rieff by equating Uint{W/2) to Uef f^lW^ / A^qK, where 
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FIG. 1: Normalized effective vortex density risff/n vs nW^ 
from Eq. ([14}. 



which is plotted in Fig. [T] Since numerical calculations 
of Uint {x) as approximated above show that 
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(15) 



we can use the argument that led to Eq. pUj) to state 
that the balance between vortex generation and escape 
occurs when 

Ba-n,ff^o = BK^im^. (16) 

The value of n corresponding to each value of Ba then 
can be obtained from Eq. or Fig. [TJ 



III. MEASUREMENTS ON YBCO STRIPS 

We performed SSM^l measurements on YBCO strips. 
A sample was prepared on a SrTiOs substrate with a 
pulsed laser deposited 200 nm thin film of YBCO. The 
sample was structured into strips varying from 2-50 /im 
in width by Ar ion etching. The SSM, in which the sam- 
ple was cooled, was placed in a liquid helium bath cryo- 
stat with three layers of ^-metal shielding. The SQUID 
used in the SSM had a pickup loop which was defined 
by focused ion beam milling and had an effective area of 
10-15 /xm^ during imaging. A magnetic induction per- 
pendicular to the sample was produced by a solenoid coil 
which was placed around the sample and SQUID. After 
the desired magnetic induction was applied the sample 
was cooled to 4.2 K and the sample was scanned. Many 
different field values have been applied to the sample dur- 
ing cooling to determine the critical induction for the 
various strip widths. The sample was warmed up to well 
above Tc between different cooling cycles. 
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a) b) 




FIG. 2: Scanning SQUID microscope images of 35/^m wide 
YBCO strips cooled in magnetic inductions of a) 5 /xT, b) 10 
HT, c) 20 fiT, and d) 50 iiT. 

In Fig. [2] SSM images are displayed of 35 /im wide 
strips for several inductions from 5 to 50 f/T. The strips in 
these images are darker than their surroundings because 
of a change in the inductance of the SQUID sensor as it 
passed over the superconducting strip. The bright dots 
are trapped vortices. As the inductions increased the 
vortex density also increased until it became difhcult to 
distinguish one vortex from the other (Fig. [2ji). In Fig. 
[5^ and Fig. [^b it is clear that at low trapped vortex 
densities the vortices tended to form one single row in 
the center of the strip where the energy is lowest. In Fig. 
[2t two parallel lines have been formed, but with some 
disorder. 



A. Critical induction vs. strip width 

The results of the measurements of the critical induc- 
tion vs. strip width are displayed in Fig. [Altogether with 
the various models. The measurements were performed 
on strips varying from 6-35 /im in width. Measurements 
on strips narrower than 6 /im were unreliable because the 
critical induction was high enough to degrade the SQUID 
operation. The critical induction for 40 and 50 fim wide 
strips was smaller than the uncertainty in the applied 
induction. 

There are two data points in Fig.[3]for each strip width: 
The upper point indicates the lowest induction at which 
vortices were observed trapped in the strip, and the lower 
point indicates the highest induction at which vortices 
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FIG. 3: Critical inductions for vortex trapping as a function 
of strip width. The squares represent -Bc+, the lowest induc- 
tions in which trapped vortices were observed, and the dots 
are Be- , the highest inductions in which trapped vortices were 
not observed. The dashed-dotted line is the metastable criti- 
cal induction Bo [Eq. ([5])], the short-dashed and long-dashed 
lines are Bl [Eq. (j4|], the absolute stability critical inductions 
calculated at a depinning temperature T^p = O.OSTc, with the 
constant a = 2/7r- or a = 1/4.- The solid line is Bk [Eq. %\. 



were not observed. This provides an upper and lower 
bound for the actual critical induction. It was apparent 
from this log- log plot that the critical induction depended 
on strip width as a power law. The best chi-square fit of 
the experimental data to the two parameter power law 
Be = a4>o/W^P yielded a = 1.96 + 0.13 - 0.15 and p = 
1.98 ± 0.03, taking a doubling of the best-fit chi-square 
as the uncertainty criterion. If we set the exponent to 
be p = 2 [Be — a$o/W^^)j a one parameter fit yielded 
a = 1.55 ± 0.27. This is to be compared with a — 1.65 
for the present model [Eq. ([5])], plotted as Bk in Fig. 
131 It should be emphasized that there were no fitting 
parameters in plotting Bk- 

Comparison of the experiment with the models of Eqs. 
([5]) and (O is straightforward, since they are only de- 
pendent on the strip width. In order to evaluate Eq. Q 
one must make an estimate of the temperature at which 
vortex freezout occurs because of the temperature depen- 
dence of ^. The depinning temperature TjTc = 0.98 used 
in Fig. 121 for both Bl curves was calculated by Maurer et 
alf^S for YBCO. In addition we used £,ybco{0) = 3 nm, 
a critical temperature of Tc = 93 K, and the two-fluid 
expression for the temperature dependence of the coher- 
ence length, resulting in £,{Tdp) — 10.39 nm. To the best 
of our knowledge the depinning temperature of YBCO 
has never been determined experimentally. Analysis of 
Eq. m shows that a Tdp closer to Tc could give better 
agreement between theory and experiment for some strip 
widths. However the difference in slopes between theory 
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and experiment becomes larger for higher T^p, making it 
appear unUkely that this is the correct model for our re- 
sults. The dependence of the Likharev model predictions 
on Tdp is displayed in Fig. 2] for a — 2/7r. For lower T/Tc 
ratios the curve moves further away from experiment. 
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FIG. 5: Comparison of experimental results on the critical 
induction for vortex exclusion in thin film niobium stripsi^ 
with various theories, labeled as in Fig. [S] 



FIG. 4: Variation of the prediction of Eq. Q (using a — 
2/7r) for the vortex exclusion critical induction on depinning 
temperature (dashed lines). The solid line is Bk [Eq. ©]. 

We also compare results of the present model with pre- 
vious work on Nb strips by Stan et alfi^ using SHPM. 
This paper reported critical inductions for 3 different 
strip widths: 1.6, 10 and 100 /im. The critical inductions 
have been compared to the various models in Fig. [5l The 
depinning temperature of T/Tc = 0.9985 used in this 
figure was experimentally determinediiS Using £,Nb{0) = 
38.9 nm results in the value ^Nb{T = Tdp) = 320 nm used 
for the Bl curves in Fig.O A reasonably good agreement 
exists between the measurements and the predictions of 
the present model. 



B. Trapped vortex density as a function of applied 
induction 

In Fig. [S] the experimentally determined density of 
trapped vortices as a function of induction for two strip 
widths is displayed. This density depends nearly lin- 
early on the difference between the induction and the 
critical induction, with a slope nearly $0 ^, in agreement 
with previous work on Nb strips by Stan et alfiS, The 
35 /im strip width data can be fit to a linear depen- 
dence of the vortex density n on Ba with a slope of 
(3.86 ±0.08) X 10i*(Tm2)-i = (0.83 ± 0.02)$o S with 
an intercept of 3.8 ± 1.3/iT. The dashed and solid lines 
in Fig. [HI are the prediction of the present model [Eqs. 
(fTTj) and (|16p respectively] without any fitting parame- 
ters. Reasonable agreement exists between the present 
model and measurements. In the case of the 6 /zm strips, 
there is an apparent saturation in the vortex density for 
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FIG. 6: Plot of the number density of vortices trapped in 
YBCO strips 35 /im and 6 fim wide as a function of mag- 
netic induction (dots). The dashed lines are the predictions 
of Eq. pip , without any adjustable parameters. The solid 
lines indicate the predictions of Eqs. (|14|l and (I16p 



inductions higher than 130 f/T. This may, however be an 
artifact due to the finite resolution of our SQUID sen- 
sor. The direction of the applied induction was reversed 
for three points in the W = 35 fim strip data to check 
for an offset in the applied induction. Such an offset, if 
present, was small, as indicated by the symmetry of the 
data around zero induction. 



7 




5 10 15 20 25 30 35 
X [urn] 



5 10 15 20 25 30 35 
X [urn] 



B=18nT 



m 

5 10 15 20 25 30 35 
X [(im] 



!_-20nT" 



^3 

5 10 15 20 25 30 35 
x[nm] 




100 200 300 400 
Ay [(im] 



60 

2; 40 

20 




B_ = 1 (iT 
Man = 12.9 nm 
Stddev = 3.29 (im 







Ay [lira] 



= 5 tiT 
Mean = 53.4 |am 
Stddev = 20.1 ^m" 



Ay [urn] 



Ij 

















10 15 20 25 



FIG. 8: (a-c) Histograms of the longitudinal spacing between 
vortices trapped in a 35 wide YBCO strip for selected 
inductions, (d) Plot of the standard deviation of the distri- 
bution of longitudinal spacings, divided by the mean of this 
distribution, as a function of induction. The relative widths of 
the distributions became narrower as the induction increases, 
indicative of ordering in a single row, until at a critical in- 
duction of about 10 there was an abrupt increase in the 
relative width as two rows started to form. 



FIG. 7: Histograms of the probability of trapping as a func- 
tion of the lateral vortex position in a 35 fim wide YBCO 
strip at various inductions. At low inductions the vortices 
trapped in a single row near the center of the film, but above 
a induction of about 10 fiT they started to reorder. At a in- 
duction of 18 jj,T the vortices were trapped in two relatively 
well defined rows. 



C. Vortex spatial distribution 

The local minimum in the Gibbs free energy at W/2 
of Eq. ^ makes it energetically favorable for vortices to 
be trapped in the center of the strip. However, as the 
vortex density increases, the vortex-vortex repulsive in- 
teraction makes it energetically more favorable to form 
an Abrikosov-like triangular pattern. Simulations on the 
trapped vortex position in strips was described by Bron- 
son et al.— In particular they predict that there should be 
a single line of vortices for inductions Be < Ba < 2A8Bc- 
Above this induction range a second line of vortices is 
predicted to form. As the induction is increased further 
additional lines of vortices are predicted to form into a 
nearly triangular lattice. 

We have investigated the distribution of vortices 
trapped in our strips at various inductions. As can be 
seen from the images of Fig. [21 even though there was sig- 
nificant disorder in the vortex trapping positions, there 



was also some apparent correlation between the vortex 
positions. An example can be seen in Fig. [71 where a 
histogram is displayed of the lateral positions of vortices 
trapped in the 35 /im wide strip for several inductions. 
At low inductions, the vortex lateral position distribution 
peaked near the center of the film because the vortices 
were aligned nearly in a single row. At a second critical 
induction of Bc2 = 11 ± 1/iT the distribution started to 
become broader. At 18 f/T there were two clear peaks in 
the distribution, corresponding to two rows. Using the 
value of Be = 3.8 ± 1.3/iT for the critical induction of 
the 35 /im wide strips from our linear fit of the vortex 
density vs. applied induction curve of Fig. [SI we found 
Bc2 = (2.89 + 1.91 - 0.93)Bc. This is consistent with 
the prediction of Bc2 = 2A8Bc of Bronson et ali^ In 
the same paper the critical induction for the transition 
from the two-row to the three-row regime is given to be 
Bc3 = 4:MBa. This gives B^s = 18.77 ± 6.42/iT using 
the same value for Be- In our measurements we saw no 
evidence for a three row regime. It was not possible to 
perform analysis at higher fields than reported here be- 
cause of limitations to the spatial resolution of the SSM. 

We also saw evidence for longitudinal ordering. In 
Fig. [5^-c histograms are displayed of the longitudinal dis- 
tances Ay between vortices in the 35 /im wide strip for 
various inductions. As expected the inter-vortex spac- 
ing distributions became narrower as the inductions in- 
creased, since the vortex mean spacings decreased. How- 
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ever the distributions became narrow faster than their 
means as the induction was increased, indicative of lon- 
gitudinal ordering, until the second critical induction _Bc2 
of approximately 10 fiT was reached. At that induction 
the relative distribution width 6{Ay) / <Ay> has a dis- 
continuous jump as a second row starts to form. A similar 
decrease in the relative longitudinal distribution width 
with increasing induction is observed in the 6 /ini wide 
strip, although the spatial resolution of the SSM was not 
sufficient to resolve vortices at the second critical induc- 
tion for this width. 

In theory there should be longitudinal ordering inde- 
pendent of the magnetic induction. After all, the Gibbs 
free energy is independent of the position along the strip 
and the only interaction that plays a role is the inter- 
action between the vortices. Differences in longitudinal 
ordering as a function of the magnetic induction could 
arise from local minima of the Gibbs free energy caused 
for example by defects in the material. For relatively low 
inductions vortices can easily be trapped in defects since 
the interaction between the vortices is small because the 
separation between the vortices is large. For higher mag- 
netic inductions the number of vortices and likewise the 
interaction between the vortices increases. This could 
mean that the vortices are more likely to trap at positions 
determined by the minimization of the vortex-vortex en- 
ergy than at positions determined by local defects. 

IV. CONCLUSIONS 

Experiments on vortex trapping in narrow YBCO 
strips using a scanning SQUID microscope, as well as pre- 



vious measurements on Nb^S, showed a critical induction 
for the onset of trapping and a dependence of the vortex 
density on the induction which were in good agreement 
with a new model which takes into account the energy re- 
quired to generate a vortex- antivortex pair. In addition, 
at low inductions the vortices formed a single row, with 
longitudinal ordering as the inductions increased. For- 
mation of a second row was observed at a second critical 
induction consistent with numerical modeling. 
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